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Quadratic fluctuations require an evaluation of ratios of 
functional determinants of second-order differential operators. 
We relate these ratios to the Green functions of the opera- 
tors for Dirichlet, periodic and antiperiodic boundary condi- 
tions on a line segment. This permits us to take advantage of 
Wronski's construction method for Green functions without 
knowledge of eigenvalues. Our final formula expresses the ra- 
tios of functional determinants in terms of an ordinary 2x2 
-determinant of a constant matrix constructed from two lin- 
early independent solutions of a the homogeneous differential 
equations associated with the second-order differential opera- 
tors. For ratios of determinants encountered in semiclassical 
fluctuations around a classical solution, the result can further 
be expressed in terms of this classical solution. 

In the presence of a zero mode, our method allows for a 
simple universal regularization of the functional determinants. 
For Dirichlet 's boundary condition, our result is equivalent to 
Gelf and- Yaglom ' s . 

Explicit formulas are given for a harmonic oscillator with 
an arbitrary time-dependent frequency. 

I. INTRODUCTION 

Evaluation of Gaussian path integrals is needed in 
many physical problems, notably in all semiclassical cal- 
culations of fluctuating systems. Typically, we are con- 
fronted with a ratio of functional determinants of second- 
order differential operators For Dirichlet boundary 
conditions encountered in quantum mechanical fluctua- 
tion problems, a general result has been found by Gelfand 
and Yaglom Working with time-sliced path integrals, 
they reduced the evaluation to a simple initial- value prob- 
lem for the homogeneous second order differential equa- 
tions associated with the above operators. The func- 
tional determinants are directly given by the value of 
the solutions at the final point. Unfortunately, Gelfand 
and Yaglom's method becomes rather complicated for the 
periodic and antiperiodic boundary conditions of quan- 
tum statistics (see Section 2.12 in and has there- 
fore rarely been used. Several papers have studied the 
functional determinants of second-order Sturm-Liouville 
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operators with periodic boundary conditions []3|- |fj, and 
related them to boundary- value problems. The calcu- 
lated determinants are all singular and were regularized 
with the help of generalized zeta- functions Q. This has 
the disadvantage of a physical quantity depending un- 
necessarily on the analyticity properties of generalized 
zeta-functions. Moreover, the auxiliary boundary-value 
problems were formulated in terms of first-order oper- 
ators, rather than the initial second-order one, making 
the treatment of a zero mode of operator with periodic 
boundary conditions unclear, and requiring additional 
work Jj|. 

In this paper we shall avoid the above drawbacks by 
developing a simple and systematic method for finding 
ratios of functional determinants of second-order differ- 
ential operators with Dirichlet, periodic and antiperi- 
odic boundary conditions. By focussing our attention 
upon ratios instead of the determinants themselves, we 
avoid the need of regularization. The main virtue of our 
method is that it takes advantage of the existence of 
Wronski's simple construction rule for Green functions. 
This permits us to reduce the functional determinants to 
an ordinary constant 2x2 -determinant formed from so- 
lutions of homogeneous differential equations associated 
with the differential operators. For semiclassical fluctu- 
ations around a classical solution, our final result will 
be expressed entirely in terms of a classical trajectory. 
Furthermore, for fluctuation operator with a zero mode, 
a case frequently encountered in many semiclassical cal- 
culations, we find a simple universal expression for the 
regularized ratio of determinants without the zero mode. 

II. BASIC RELATIONS 

The typical fluctuation action arising in semiclassical 
approximations has the form 

l-tb rtt, j[,r 

A[x] = / dtL(x(t),x{t)) = / dt— [x 2 - n 2 (t)x 2 ] . 
Jt a Jt a 2 

(1) 

The time-dependent frequency £l(t) can be expressed in 
terms of the potential V(x) of the system as 

Q 2 (t) = V"(x cl (t))/M, (2) 

where x c \(t) is a classical trajectory solving the equation 
of motion (for examples see Section 17.3 of Ref. 0). 
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Mx 



-V\x). 



(3) 



The action (P describes a harmonic oscillator with a 
time-dependent frequency Q(t). For this system, both 
the quantum mechanical propagator and the thermal par- 
tition function contain a phase factor exp[i^4 c i], where 
A c \ = A[x c \] is the action of the classical path x c \(t). 
The phase factor is multiplied by a fluctuation factor 
proportional to 



F{t b ,t a ) 



/DctKi 
V Dct/f 



-1/2 



(4) 



where K\ = —df — f2 2 (i) = K — Q 2 (t) is obtained as 
the operator governing the second variation of the action 
A[x] along the classical path x c i(t): 



S 2 A[x cl ] 
5x(t)5x{t r 



= 8{t-t')K 1 . 



(•5) 



The ratio of determinants (|l|) arises naturally from the 
normalization of the path integral Q and is well-defined. 
The linear operator K\ acts on the space of twice differ- 
cntiablc functions y(t) = 5x(t) on an interval t G [t a ,tb] 
with appropriate boundary conditions. In the quantum- 
mechanical fluctuation problem, these are Dirichlet-like 
with y(tf,) = y(t a ) = 0, in the quantum-statistical case, 
they are periodic or antiperiodic with y(tb) — ±y(t a ) and 
y(tb) = iy(iet)- The operator K in the denominator of 
(|j) may be chosen as K$ (Dirichlet case) or Kf = Kq—uj 2 
(periodic and antiperiodic cases), respectively, where to 
is a time-independent oscillator frequency. Then the op- 
erator K is invertible, having the Fredholm property 



DetiCi 
DetK 



(6) 



(a possible multiplicative anomaly being equal to unity 
Furthermore, since the operator K~ X K\ is of the 
form I + B, with B an operator of the trace class, it has 
a well-defined determinant without any rcgularization. 

To calculate F(tf,,t a ), we introduce a one-parameter 
family of operators 



Kg = 



(7) 

depending linearly on the parameter g G [0,1], and re- 
ducing to the initial operator K\ for g = 1. Then we 
consider the eigen- value problem 



K g(t)y n (g;t) = X n {g)y n (g;t), 



(8) 



with eigenvalues X n (g)- The eigenfunctions y n (g;t) sat- 
isfy the orthonormality and completeness relations 



dty n (g;t)y m (g;t) = 5 nm , 
^2 yn ^ 9 ' *') = s ( l - t ')- 



(9) 
(10) 



The completeness relation permits us to write down im- 
mediately a spectral representation for the Green func- 
tion G g (t, t') associated with the differential equation (|§). 
By applying K g (t) to 



G S (M') = E 



yn(g;t)y n (g;t') 



(11) 



and using fl|) , ( |To| ) , we verify the validity of the defining 
differential equation 



K g (t)G g (t,t') = 5(t-t'). 



(12) 



In terms of the eigenvalues X n (g), the determinant 
would read 



Dct K-^g = C Yl 



A»(g) 



(13) 



where C = Det (K ~ 1 Kq) is a constant of the g- 
integration, which still may depend on tb,t a . Since the 
infinite product of ratio of the eigenvalues X n in Eq. ( p^ ) 
converges uniformly for all g G [0, 1], we can differentiate 
this equation to obtain 



dglogBctK^Kg = 



KM 



(14) 



Differentiating Eq. (ph, and using the condition ( |10[ ) gives 
for all boundary conditions, 



X' n (g) = - dtn 2 (t)y 2 n (g;t). 



(15) 



This may be inserted into ([l4]). Because of the conver- 
gence of sum in (|l4|), summation and integration can be 
interchanged, and using the spectral representation ( p"l| ) 
we find the compact formula 

d g log Det (K- l K g ) = -Tr [n 2 {t)G g (t, t')] . 



dtn 2 (t)G g (t,t) (16) 



By integrating this equation in g, we obtain the ratio of 
functional determinants in the form 



Deti^T 1 !^ = Ccxp<^ - / dg' l dtVL 2 {t)G g ,{t,t) 



(17) 



with the same integration constant C as in Eq. (|13|). It 
is fixed calculating the same expression for 5 = where 
the left-hand side is well-known. In the case of Dirichlet 
boundary conditions where K = Kg, the left-hand side 
is trivially unity. For periodic and antiperiodic boundary 
conditions where we take K = Kq — lo 2 = —d 2 — uj 2 , the 
most convenient way to normalize the right-hand side is 
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to go to g — 1 and choose the frequency fi 2 (t) to be equal 
to the constant frequency lo 2 . The left-hand side is again 
unity thus fixing C. 

Having determined C we set g = 1 in Eq. ( |l7|) and 
obtain the final result for the operator K\. In the sequel 
we shall evaluate the right-hand side of formula ((lj) using 
explicit Wronski constructions of the Green function for 
the different boundary conditions. 



III. WRONSKI'S CONSTRUCTION OF GREEN 
FUNCTIONS 

The general solution of the differential equation ( |l2] ) 
may be expressed in terms of retarded and advanced 
Green functions which have the general form 



G g -\t,t') = G g + \t',t) = e tt > A g (t,t'), 



(18) 



where ®tt' = — t') is Heaviside's step function which 
vanishes for t < t' and is equal to unity for t > t' . The 
function A g (t, t') satisfies the homogeneous differential 
equation corresponding to ( |l2| ) . This is seen by applying 
the operator K g {t) to ( |18] ) and making use of the identity 
tS'(t) = -S(t): 



K g (t)G g -\t,t')=e tt ,K g (t)A g (t,t') 



(*-*') 



v g\ v )^g 
-2d t A g (t,t') 



S(t-t'). 



(19) 



Since the right-hand side must be equal to 5(t — t'), the 
function A g (t,t') has to satisfy the homogeneous differ- 
ential equation 

K g (t)A g (t,f) = 0, for t>t' , (20) 

while the bracket in ( |l9|) must be equal to 1 at t = t'. 
Upon expanding A g (t, t') around t = t' , this leads to the 
conditions 

A g (t,t)=0, d t A g (t,t')\ t , =t = -I. (21) 

Equation ( p0| ) is solved by a linear combination 

A g (t,t')=a g (ti) Vg (t)+p g (t')t g (t) (22) 

of any two independent solutions ri g (t) and £, g {t) of the 
homogeneous equation 



K g (t)h g (t) = [~d 2 - g n 2 (t)] h g (t)=Q. 



(23) 



Their time-independent Wronski determinant W g = 

Vg£,g ~ Vg£g is nonzero, so that we can determine the 
coefficients in the linear combination (]2^) from ( pl| ) and 
find 



A g (t,t') 



iv 9 m g (t')-ut)v g (t')] = -A g (t',t). 



The right-hand side contains the so-called Jacobi com- 
mutator of the two functions r) g (t) and £ g (t). Here we 
list a few algebraic properties of A g (t, t') which will be 
useful in the sequel: 

a f . _ A g (t b ,t)A g (t',t a ) - A g (t b ,t>)A g (t,t a ) 



Ag(t b ,t)d tb A g (t b ,t a ) ~ Ag(t,t a ) = Ag (t fe , t „ ) d t Ag (t fe , t ) , 

(26) 



A a (i, i )ft„ A 9 (4, *a) + t) = A g (t b , t a )d t A g (t, t a ). 

(27) 

Note that the solution ([l^) is so far not unique, leaving 
room for an additional general solution of the homoge- 
neous equation (|23j) 



G g (t, = 9«, A g (t, t') + a,(f )»&,(*) + &,(*%(*) (28) 

with arbitrary coefficients a g (t') and b g (t'). This ambi- 
guity is removed by appropriate boundary conditions. 

Consider first the quantum mechanical fluctuating 
problem with Dirichlet boundary conditions y{g]t b ) = 
Ug(ta) = for the eigenfunctions y{g\t) of K g , implying 
for the Green function the boundary conditions 



Gg(t b ,t)=0, t b ^t, 
Gg(t,t a )=0, t^t a . 



Substituting ( |28| ) into 
pair of equations 



(29) 

|) leads to a simple algebraic 



a g (t)Vg(t a )+b g m g (Q=0, (30) 
a g (t)Vg(t b ) + b g (t)£ g {t b ) = -A(t b ,t). (31) 

We now define a fundamental matrix A g as the constant 
2 x 2-matrix 



A | £ff(*a) 

9 \%(t b ) Uh) 



and observe that under the condition 

&etkg = WgAg{t a ,t b )^<d, 



(32) 



(33) 



the system (|3l]) has a unique solution, so that the co- 
efficients a g (i) and b g (t) in the Green function ( p8| ) are 
easily calculated. Making use of identity (|25|), we obtain 
Wronski's well-known formula 



Gg{t,t'): 



e tt >Ag(t b ,t)A g (t',t a )+Q tlt A g (t b ,t')Ag(t,t a ) 



Ag(t a ,t b ) 



(34) 



(24) 



For Dirichlet boundary conditions, this equation yields 
a unique and well-defined Green function assuming the 
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absence of a zero mode of the operator K\ with these 
boundary conditions. Such a mode would cause prob- 
lems since r]i(t a ) = T)x(t b ) = would make detAi = 0, 
thus destroying the property ([33j) which was necessary to 
find ([mJ) . Indeed, the Wronski expression (124) is undeter- 
mined since the boundary condition r?i(t J = together 
with (|30|) imply fi(i ) = 0, making Wi = r)i£j. - r)x(i 
vanish at the initial time t a and thus identically in t. 

Consider now the quantum statistical fluctuation prob- 
lem with periodic or antiperiodic boundary conditions 
y{g;t b ) = ±y(g;t a ), y(g;t b ) = ±y(g;t a ) for the eigen- 
functions y(g;t) of the operator K g {t). For the Green 

function G g (t,t'), these imply 



Gg (tb,t') — ±Gg (t a ,t), 

P P 
Gg(t b ,t')=±Gg(t a ,t'). 



(35) 



In both cases, the frequency Q(t) and the Dirac delta 
function in Eq. (jl^) are also assumed to be periodic or 
antiperiodic in time with the same period. Inserting ( ps| ) 
into ( |35| ) gives now the equations 

a(*)(r? 6 TVa) + &(*)(& T £«) = ~A(t b , t), 

o(t)(ifc T Va) + 6(t)(6 T L) = -d t A(t b , t). (36) 

For brevity, we have omitted the subscripts g and written 
£a,b,ria,b for £g{t a ,b),'ng{'ta,b)- Defining now the constant 
2x2 -matrices 



Vb TVa 6> T £o 
Vb TVa ib T ia, 



(37) 



IV. MAIN RESULTS AND RELATION TO 
GELFAND-YAGLOM'S INITIAL- VALUE 
PROBLEM 

Excluding at first zero modes, we evaluate formula (|1^ 
for ratios of functional determinants. The temporal inte- 
gral on the right-hand side can be performed efficiently 
following Ref. @. Here we present an even more di- 
rect method, by which we express the result in terms of 
solutions of Gelfand-Yaglom's initial-value problem for 
Dirichlet boundary conditions, and of a dual problem for 
periodic and antiperiodic boundary conditions. 

Dirichlet Case 

The Gelfand-Yaglom initial- value problem consists in the 
search for a function D g (t) solving the following equa- 
tions: 

K g {t)D g (t)=0; D g (t a ) = 0, b g (t a ) = l. (41) 

By differentiating these three equations with respect to 
the parameter g, we obtain for D' g {t) = d g D g (t) the in- 
homogeneous initial-value problem 



K g (t)D' g (t) = n 2 (t)D g (t)- D' q (t a ) = 0, D' q (t a ) = 0. 



(42) 



the condition analogous to (p3|) 



The unique solution of equations (|4l]) can easily be ex- 
pressed in terms of our arbitrary set of solutions rj g {t) 
and £ g (t) as follows 

D g (t) = ^m-^aW) = Ag{ ^ t) (43) 
W a 



detAg = Wghg (t a ,t b ) ^ 



(38) thus leading to 



with 



A g (t a , t b ) = 2± d t A g (t a , t h ) ± d t A g (t b , t a ) 



(39) 



enables us to obtain the unique solution to Eqs. (36) 
After some algebra using the identities 
the expression (E8h for Green functions wit 



and (27) 



antiperiodic boundary conditions 
the form 



<(M') = G 9 (M')T 

[A g (Ma) ± A g (t b ,t)][A g (t',t a ) ± A g (t b ,t')} 

p 

Ag(t a ,t b )A g (t a ,t b ) 



i periodic and 
5|) can be cast into 



(40) 



The right-hand side is well-defined unless the operator 

K\ has a zero mode with rj b = irj a , Vb = ^Va, which 

p 

would make the determinant of the 2x2 -matrix A g 
vanish. 

Note that the Green functions (Q) and (40|) are both 
continuous at t = t', as is necessary for calculating the 
associated ratios of functional determinants from formula 
@, which we shall now do. 



D g it b ) - 



Wg ^ 



= A g (t a ,t b ). 



(44) 



In terms of the same functions, the general solution of the 
inhomogeneous initial-value problem ([42| ) can be seen to 
have the form 



D' g (t) = I <ft / n 9 (t / )A fl (i,t , )A fl (t 01 i / )- (45) 



Comparison with (p4|) shows that at the final point t = t b 



D' g (t b ) = -A g (t a ,t b ) J dtn 2 {t)G g (t,t). (46) 

which together with (E3) implies the following simple re- 
lation for the Green function (|34] ) with Dirichlet 's bound- 
ary conditions: 

Tr[f7 2 (i)G 5 (t,i')] = -d g \og = -dg\ogDg(t b ). 

(47) 
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Inserting this into (|T7|) , we find for the ratio of functional 
determinants the simple formula 



(48) 



VebKo l K B = CD g (t b ). 



The constant of integration is fixed by applying ( p8[ ) to 
the trivial case 5 = 0, where Kq — ~df and the solution 



to the initial- value problem (41) is 
D (t)=t-t a . 



(49) 



At g = 0, the left-hand side of ( J50| ) is unity, determining 
C = (tb — ta)^ 1 and the final result for g = 1: 



detAi /DetA £>i(4) 



Wo 



th — t„ 



(50) 



This compact formula was first derived by Gelfand and 
Yaglom g via a direct calculation of the determinant 
arising in a time-sliced path integrals pL 



Periodic and Antiperiodic Case 

Our technique makes it straight-forward to derive an 
equally compact formula for periodic and antiperiodic 
boundary conditions. For this purpose we introduce an- 
other homogeneous initial-value problem whose bound- 
ary conditions are dual to Gelfand and Yaglom's in (41): 



K g {t)D g (t) = 0; D g (t a ) = l, D g (t*) = 0. (51) 

In terms of the previous arbitrary set T] g (t) and t; g (t) of 
solutions of the homogeneous differential equation, the 
unique solution of (|5l|) reads 



D g (t) 



(52) 



This can be combined with the time derivative of (^) at 
t = t b to yield 



Dg(t b ) + Dg(tb) =i[2- Ag (t«,t 6 )]. 



(53) 



By differentiating Eqs. ( |51| ) with respect to g, we ob- 
tain the following inhomogeneous initial-value problem 
for D' g (t) = d g D g {t): 

Kg(t)D' g (t) = n 2 (t)D' g (t); D' g (t a ) = 0, D' g (t a ) = 0, 

(54) 



whose general solution reads in analogy to 



D' g (t) = -J dt'n 2 (t')A g (t,t')A g (t a ,t'), (55) 

where the dot denotes the time derivative with respect 
of the first argument of A g (t,t'). With the help of iden- 
tities (||) and @, the combination D'(t) + D' g (t) at 
t = tb can now be expressed in terms of the periodic and 
antiperiodic Green functions (E3) , in analogy to ([46]) , 



D'(t b )+D'(t b ) = ±Ag(t a ,t b ) / dtn 2 (t)G g (t,t) 



(56) 

Together with (|53|), this yields for the temporal integral 
on the right-hand sides of (|l6|) and (|l7|) the simple ex- 
pression analogous to (E?j) 



Tr[n 2 {t)G g {t,t') 



., , , detAg 

-dg fog ' 



-9g fog 



Wg 

2TD g (tb)TDg(t b ) . (57) 



This is inserted into formula Eq. (|17|) yielding for periodic 
and antiperiodic boundary conditions 



C 



2TDg{tb)TD g {t b ) 



(58) 



where K = Kq — uj 2 = — d 2 — u) 2 . The constant of in- 
tegration C is fixed in the way described after Eq. (|rj|). 



We go to g — 1 and set fr(i) 



For the operator 



Kf = —d 2 —u> 2 , we can easily solve the Gelfand- Yaglom 
initial- value problem (^Tj) as well as the dual one (|l]) by 

Df(t) = - wa[u{t - t a )], D<?(t) = cos[u;(t - Q], (59) 



so that (^8|) determines C by 



1 = C 



4sin 2 [w(t(, — t a )/2] periodic case, 
4cos 2 [cj(ih — t a )/2] antiperiodic case. 



(60) 



Hence we obtain the final results for periodic boundary 
conditions 



Det( A-.A- l) = ^/^ 



= 2 - D^tb) - Diitb) 
" 4sin 2 [a;(t 6 -t a )/2] : 

and for antiperiodic boundary conditions 

detA? /DetA^ a 



(61) 



Det (K^Kx) = 



Wi I W? 

2 + b^h) + Dijtb) 
4cos 2 [oj(t b -t a )/2} ' 



(62) 



The intermediate expressions in ©, ©, and (H 
show that the ratios of functional determinants are ordi- 
nary determinants of two arbitrary independent solutions 
rji(t) and £i(t) of the homogeneous differential equation 
K iy (t) = [-Of - fl 2 (t)]y(t) = 0. As such, the results 
are manifestly invariant under arbitrary linear transfor- 
mations of these functions (?7i,£i) — ► (?7i,£i)- 
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V. EXPRESSIONS IN TERMS OF CLASSICAL 
TRAJECTORY 

In semiclassical fluctuation problems, the time- 
dependent frequency fl 2 (t) is determined by the clas- 
sical solution x c i(t) of the equation of motion (||) via 
Eq. (g). In this case, the above results can be made 
quite explicit by expressing the solutions D\(t) and D\(t) 
of the initial-value problems (^) and ( pl| ) directly in 
terms the classical trajectory x c \(t) if this is specified 
in terms of its initial position x a and initial velocity x a 
as x c \(t,x a ,x a ). Given such a trajectory x c \{t, x a , x a ) = 
x a £>i(t) + x a Di(t) the solutions of ([4l| ) and J5l| ) can be 
written in the form 



D 1 (t) 



dx c \(t) X a7 X a ) 



Di{t) 



3x c \(t : X a: X a ) 
dx a 



(63) 



As an example, take a harmonic oscillator where for- 
mulas (|63[) are given explicitly by the previous expres- 
sions (|59]). For a classical path, we can use the equation 
of motion ([}]) and a partial integration to express the 
action as a surface term 



A[x c \\ = M(x b x b - x a x a )/2, 



where 



Xb 

x b 



X a Di(t b ) + XaD^tb), 
X a D 1 (t b ) + X a Dl(t b ). 



(64) 



(65) 



With the help of Eqs. (65), we can write the action ( |64| ) 
as a function of initial and final positions x a and x b , and 
of the time difference t b — t a : 



Ac\ {x a , X b: t b t a ) 

x[bi{%)xl 



M 



2Dx{t b ) 

2x b x a + Di(t b )x 2 a ] 



From this we obtain directly 

d 2 Aci{x ai x bl t b 



£>i(i 6 ) = -M 



ta) 



dx a dx b 



(66) 



(67) 



so that the ratio (|50|) of functional determinants for 
Dirichlet boundary conditions becomes 

DetK' 1 ^ = -M 



d 2 Aci(x a ,Xb,t b - t a ) 



dx a dx b 



(h-t a ). 
(68) 



The right-hand side is known as one-dimensional Van 
Vleck-Pauli-Morette determinant (see Section 4.3 in 0). 

In the case of periodic and antiperiodic boundary con- 
ditions, we find from Eq. (pq) 



(69) 



2TD 1 {t b )TD 1 (t b ) = 2± 

d 2 Acl(x a ,X b ,tb - t a ) d 2 Acl(x a ,Xb,tb - t a ) 



d 2 A c i(x a ,x b ,t b - t a ) 



dx a dx b 



dx 2 



dx 2 



which determines the ratio of functional determinants 
(|S8|) in terms of the classical action, in analogy to (|6§|). 
For a harmonic oscillator with the classical action 



Acl(x a: x bl t b t a ) 



Muj 



2 sinui(t b — t a ) 



x [{x 2 b + x 2 a ) cos Lo(t b - t a ) - 2x b x a ], (70) 
and we obtain Di(tb) = ui~ x sinuj(tb — t a ) as in (|59|) and 



U(f n fj u J(f ] _Usm 2 [uj(tb-t a )/2} , , . 
in agreement with the previous results (|50|), (^TJ), and 



VI. TREATMENT OF ZERO MODE 

Consider now the often encountered situations that the 
operator K\ has a zero mode. In path integrals, such 
a zero mode arises for example from the translational 
invariance along the time axis of a classical solution in 
a potential V(x). As in the last section, the squared 
frequency £l 2 (t) is determined by (||). 

For simplicity, we shall assume the presence of only a 
single zero mode, which we choose as one of two indepen- 
dent solutions of the homogeneous differential equation, 
say rj{t). For Dirichlet boundary conditions, we call this 
a Dirichlet zero mode, satisfying 



Vb = 0, i] a = 0. 



(72) 



For periodic and antiperiodic boundary conditions, the 
zero mode satisfies 



Vb T Va = 0, f] b T Va = 0, 



(73) 



respectively. As pointed out earlier, the Wronski con- 
struction for evaluating ratios of functional determinants 
is not applicable here since the conditions j33|) and d3§| ) 
are violated as a consequence of ([72]) and (|73|) . In or- 
der to enforce ([IJ]) and (|38|), we modify the boundary 
conditions for eigenfunctions y (t) of the operator K\ by 
a small regulator parameter e > 0, and determine new 
eigenfunctions y e (t) with y € {t) — » y(t) and A e — » A for 
e — > 0. The specific form of regularized boundary condi- 
tions will be irrelevant. It is merely required to keep the 
boundary-value problem self-conjugated. For instance, 
the Dirichlet boundary may be slightly modified to 

v l-er,l = 0, r]l + ef]l = Q, (74) 

the periodic and antiperiodic ones to 

tjI — ± cosh er]l + sinh e rjl, 

rf a = sinhe r\\ ± coshe rj\ . (75) 



Whereas the zero mode ?y (t) satisfies (|72j) or ( |73| ) , the 
modified function ^(t) is no longer a zero mode, but has 
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an eigenvalue SX e of K\, which goes to zero for e — > 0. As 
long as e is nonzero, the Wronski construction provides 
us with a regularized determinant Det K\ which tends 
to zero in the limit e — > 0. In terms of the independent 
solutions rj{t) and of K\y = 0, this determinant is 
given for the regularized Dirichlet boundary conditions 
(f74|), to first order in e, by 

DetK{ = Det Ki + ^ (r?„6 - »7a& + VbL - Vb£a)- (76) 

The determinant Deti^i vanishes, and the constant 
Wronskian 



W = rja^a - fl a ^ a = T] b £ b - l) b £ b 



is, by (72), equal to 
W = 



(77) 



(78) 



Simplifying (T7q) further with the help of (172), we obtain 



V>etK{ 



W 



w 



rUbivl + vl)- (79) 



For the regularized periodic and antiperiodic boundary 
conditions (^q) , the determinant reads, to first order in 
e: 

IkStKl = ^ (rjbta ~ Vatb - Vbia + Va£b) , (80) 

with the same Wronski determinant (^7|) whose con- 
stancy implies, together with ([73]), that 



VbiibTia) -Vbi&TL) = 0. 
Using (f73|) once more in (|80|), we find 



BetKl = =F 



W-qi 



(% 2 -% 2 )(6TCa). 



(81) 



(82) 



In order to find a finite expression for the functional de- 
terminant we must divide out the eigenvalue 5X e before 
taking e — > 0. From the regularized eigenvalue equation 

Kxrfit) = 5\ e rf(t), (83) 

with rj e (t) normalized as in (|^) , we find to first order in e 



TjKirf — (rffj — f] e rj) 



dtrf{t) = SX € 



Taking into account the regularized boundary conditions 
(||) and @ for rf(t), as well as the conditions @ 
and (|73|) for T](t), gives for the eigenvalue of the Dirichlet 
would-be zero mode r] e (t) 

<5A e = iW b - rjaVt = ~ e (W6 + VaVt), (84) 
and for periodic (antiperiodic) boundary conditions: 

SX e = fi b (ril T Vl) ~ Vb(r]b =F vl) = Te(r)bVb ~ VbVb)- ( 85 ) 



These equations enable us to remove SX e from the regu- 
larized determinants (|7^) and (jS2|). Defining the deter- 
minant without the zero mode by (see Section 17.5 in 

i). 



Detifi 

Det'Ki = lim - 1 . 

e-0 5X e 



(86) 



we obtain from (|79|) and (|84|) for Dirichlet boundary con- 
ditions 



Det'if 1 = -^%lim-^ +?) ' 



W 2 T"6 r) b ril + r) a r)e W 2 i) b f) a ' 



For periodic and antiperiodic boundary conditions, the 
result is from © and ©: 



Det'ifi 



lim 



??bVF e^o %r)| - 77 & 7?g ?7f,VF 
which by (|8l]) becomes 



Det'ifi = - 



6 T £a 



, (88) 



(89) 



Formulas ( |S7j ) and ( |39| ) are useful for semiclassical cal- 
culations of path integrals whose equations of motion 
possess nontrivial classical solutions such a s olito ns or 
instantons [EJ, as will be illustrated in Section |Vnj| . 

Note that our final expressions (p7| ) and (|89|) for the 
functional determinants are independent of the specific 
choice of regularization. 



VII. TIME-DEPENDENT HARMONIC 
OSCILLATOR 

To illustrate the power of the formulas derived in this 
work consider the time-dependent harmonic oscillator de- 
scribed by the Lagrangian (|l|). The path integral for- 
malism for such a system with the Dirichlet boundary 
conditions was studied in several papers Jl(|- fl2|| . Here 
we rederive their results and generalize them to periodic 
and antiperiodic boundary conditions. Due to the ab- 
sence of time-translational invariance of the Lagrangian 
(pj), a zero mode can be excluded here. For the Wron- 
ski construction, we take two independent solutions of 
Eq. (H) as follows 



T)(t) = q(t) COS 4>{t) , f (t ) = q{t) Sin <j>(t) 



(90) 



with a constant Wronski determinant W. The solutions 
•q(t) and £(i) are parametrized by two functions q(t) and 
q(t) satisfying the constraint 



j>(t)q 2 (t) = W. 



(91) 



The function q(t) is a soliton of the Ermankov-Pinney 
equation 13 
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q + n 2 {t)q 



W 2 c 



0. 



(92) 



For Dirichlet boundary conditions we insert (|90| ) into 
(|50|), and obtain the ratio of fluctuation determinants 
in the form 



Deti^lfi 



1 q(t a )q(t b )sm[<l>{t b ) - (/>(t a )] 



W 



tb — t a 



(93) 



For periodic or antipcriodic boundary conditions and 
Q(t), the functions q(t) and <p(t) in Eq. ( |90| ) do not in gen- 
eral have the same periodicity. This is possible because 
of the nonlinearity of Eqs. (^1J) and (|9^). Moreover, since 
we are assume here the absence of a zero mode with such 
boundary conditions, it is a necessary property of the so- 
lutions of the homogenous equations (^3|). Substituting 
© into © and (f|), we obtain the ratios of functional 
determinants for periodic boundary conditions 



-x 2 {t) = V{x{t))+E 1 



(96) 



where x (r) = dx (r) /dr, and E is the integration con- 
stant corresponding to the particle energy in the inverted 
double-well. For the splitting between ground state and 
first excited state, we must study the path integral for the 
evolution amplitude over a large but finite time interval 
( T aj T 6)- I n a semiclassical approximation, this is domi- 
nated by periodic solutions of with energy E < 0, whose 
turning points lie close to the minima of the double- well. 
We consider first a single sweep across the central bar- 
rier from a turning point at x(r a ) — x a to x(rb) = Xb, 
where the velocities vanish: x(r a ) — ifo) — 0, so that 
the energy is given by 



E = -V (x b ) = -V (x a ) = -— (x 2 - a 2 )' 



(97) 



Det^- 1 A- 1 =4sin- 2 U;fa " ta) 



x <^ 4 sin' 



W b ) - </>(t a )] 



[q(tb)q(tg) - q(t a )q(t b )} 
sm[(j){tb) - <p{t a )\ 



[q(t b ) - q(ta)} 2 
q{t a )q{t b ) 



COs[(j)(tb) - <p(t a )} 



(94) 



For antiperiodic ones, we must interchange sin — > — cos. 
By a linear combination of the solutions (90) we can al- 
ways redefine 4>{t) such that <p{t a ) = 0. 

In the literature, only formula Jj93) for the Dirichlet 
case appears to be known (see |10| - [[i"2"||). Formulas 
( |94| ) for periodic and antiperiodic boundary conditions 
are new, except for predecessors in a time-sliced formu- 
lation (see Section 2.12 in |lj). The present derivation 
is, however, much simpler than that of the predecessor 
since we have been able to take full advantage of Wron- 
ski's simple construction method for Green functions. 



For a single sweep this implies 

Xb = -x a , x b < a. 



(98) 



For an infinite time interval (t ai t^,) , the sweep connects 
the potential minima with each other, in which case 
x a = —Xb = a and E = Then Eq. ( |96| ) can easily 
be integrated yielding the well-known kink solution |lj] 
centered around some finite tq = (r a + t&)/2: 



a^ci (t) = a tanh 



w (t - T ) 



(99) 



With the explicit energy (97), the equation of motion 
(BQ) reads 



i 2 (r)^^(x 2 -x 2 ) (b 2 -x 2 ) 



(100) 



where b 2 = 2a 2 - x\ and x 2 < x\ < b 2 . Integrating 



Eq. (hOfj) gives 



VIII. FLUCTUATION DETERMINANT OF 
INSTANTON 



dt 



c(r) 



-—(r-r b ). (101) 
la 



As an application of our formulas we derive the func- 
tional determinant of the quadratic fluctuations around 
an instanton which governs the energy level splitting of a 
quantum mechanical point particle in a double- well. Set- 
ting the mass equal to unity, for simplicity, we consider 
a potential of the form (jj 



I'M -£,(-• 



(95) 



The tunneling through the central barrier is controlled 
by the solution of the equation of motion at imaginary 
time r = —it, which can be integrated once to yield the 
energy conservation law 



It is useful to introduce a normalized coordinate y(r ) = 
x (r) jxb moving between —1 and 1, and rewrite (101) as 



vi T ) 

If dt to , K 

bj V(l-^)(l-^) = ^ (T - T&)+ 6- ( 102 > 

The parameter m is equal to x\ / b 2 < 1 and determines 
the constant k on the right-hand side via the complete 
elliptic integral of the first kind 



k = K (m) 



dt 



V^l -t 2 ) (1-mt 2 ) 



(103) 
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This constant fixes the period T via formula (101) for 
t = T a as follows 



where e = E(m) is given by the complete elliptic integral 
of the second kind 



2k = —bT. 
2a 



The general solution of Eq. ( |102|) is 

x c i (t, n,m) = x b sn (z (r) ; m) , 



(104) 



(105) 



where z (r) = w6 (r — T b ) / 2a + k, so that z b = k, z a — 
—k and sn (z; m) is the elliptic function running from —1 
to 1 for t G (T a ,Tb), thus ensuring the correct boundary 
conditions x cX (r 6 ) = x fc , x cl (r a ) = -x b . 

According to Eqs. (0) a rid (p| ), the fluctuations Sx (t) = 
y (r) around the solution (llOq ) are governed by the differ- 
ential operator K\ (r) = — d 2 / dr 2 +uj 2 (3x^ — a 2 ) / 2a 2 . 
The boundary conditions are eventually irrelevant for the 
level splitting in the ground state, since this will require 
taking the limit of an infinite time interval. As an exam- 
ple, we consider here Dirichlet boundary conditions for 
eigenfunctions of the operator K\ (r) : y b = y a = 0, The 
derivation of fluctuation determinant requires, in general, 
two independent solutions of the homogeneous differen- 
tial equation which after going over from the time r to 
the variable z(t) takes the Lame's form 



h(z) 



2a 2 2 , 

— 6m sn (z; m) 

cr 



h(z) = 0. 



(106) 



By translational invariance, the first independent 
solution r\ (t) to this e quat ion is the derivative 
(d/dr b ) x c i(z(t); to) of Eq. (105). Normalizing, we have 
explicitly 



t]{t) =N 



dx c i (r, n, m) 



dr b 



. to 



= —N—bxb cn (z; m) dn (z; m) , 
2a 



(107) 



whose time derivative is 



V (t) = N [^ b ) x b sn(z;m) 

x [dn 2 (z; m) + m cn 2 (z; m)] 



(108) 



Here cn(z;m) and dn(z;m) are the elliptic functions. 
The normalization factor N is determined by the condi- 
tion (|9|) as follows 



N- 2 = xi ( 



\2a 



dz cn 2 z dn 2 z . 



(109) 



Performing the integral yields 

*r 2 4a 2 
N 



3(m + l) \2a 
x [(1 — m) k — (m + 1) e] 



(110) 



E (m) = / dl 





1 — mt 2 
1-t 2 ' 



(111) 



The solution (107) is a zero mode of the operator 
K\ (r) , since it satisfies Dirichlet boundary condition 



fib = Va 



-N — bx b cn k dn n — 0. 
2a 



(112) 



Because of the property (112), the Dirichlet determinant 
( pp| ) vanishes, so that we may only calculate the primed 
determinant according to formula (87). This does not 
require th e sec ond independent solution £ (r) of Eq.(106). 
From Eq.(108) we observe 



Vb= - 1 la= N ^—b) x b (1 - m) 



Inserting ( JL13| ) with (|11C|) into Eq. 
ately 



(113) 
we find immedi- 



Det'ATi 



1 



4a 2 [(to + 1) e - (1 - to)k] 
3a; 2 (to+1)(1-to) 2 



(114) 



Let us turn now to the limit of an infinite time interval 
where E — > and x b and 6 go to the constant a, the 
parameter to tends to unity as (1 — m) — > 16 exp (—2k), 
and e — > 1. Using Eq. (104), we obtain from (114) 



Det'JsTi 



24w 3 24w 3 



In the same limit 



sn 



2 a 



b (r - r 6 ) + fc; 1 



— > tanh 


"w(t-t )" 




2 



(115) 



(116) 



so that (105) reduces to the limiting kin k so lution ( |9 
for which the fluctuation determinant is (115). The pres 
ence of exponentially divergent factor e" T like in Eq. (115) 
is a special future of fluctuation determinants in the imit 
T — > 00. It also appears if one derives the determinant of 



harmonic differential operator K = K\ 



-d 2 /dr 2 



which governs the fluctuations around the trivial con- 
stant classical solution x c \ (r) = a, with the same Dirich- 
let boundary conditions. Indeed, inserting of the inde- 
pendent solutions r\ (r) = coshwr, £ (r) = sinhwr into 
Eq. ( |50| ) yields directly 



Det K% 



sinhwT 



2to 



(117) 



where the right-hand side being the large-T limit. Cer- 



tainl y, wh en considering the more relevant ratio of (115) 
and (117), we obtain the finite result 
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Det'if i / Det K% 



I2u 2 



(118) 



which agrees, of course, wit h pr evious calculation |lj]. 

The primed determinant ( |llq ) can also be derived from 
( p7| ) using only the asymptotic behavior of the indepen- 
dent solutions 77 (r) and £ (r) at T — > 00 |U^ |l). For 
this purpose, we set the particle energy in (|97|) equal to 
zero. The the ellipti c fun ctions degenerate into hyper- 
bolic, simplifying Eq.( 106 ) to 



h (z) - 2 (2 - 3 cosh -2 z) h (z) = 0, 



(119) 



where z (r) = to (r — to) /2 with to = (rf, + r a ) /2. Now 
we are looking for the asymptotics of two independent so- 
lutions to this equation at t\> — r a — > 00. For the solution 
corresponding to (107) we find 



r?(r) 



, T aw 
-TV — cosh 
2 



-2awA^e _2|z| , (120) 



the proper normalization factor being N~ 2 = 2a 2 u)/3. 
For the second independent solution, the asymptotic be- 
havior may be deduced from the constancy of Wronskian 
W(r},£) as follows 
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£(r) 



±e 



2\z\ 



(121) 



Here the normalization is irrelevant since the expres sion s 
© is independent of it. The solutions (|l20|) and (|j~2l]) 
have the Wronskian W(r], £) = — 4au> 2 N, and the asymp- 
totic boundary conditions for r\ (t) 



rib = r, a « -2aojN e~ w «>/ 2 , fib 
and for £ (r) 

j(r b -T a )/2 



-r) a w ujr/b, (122) 



6> — -£a 



(123) 



Inserting this into formula (g7|) , with the right-hand side 
r ewr itten as — £, a £,b/W 2 , we obtain once again the result 
(11). 
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